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Abstract 1. We consider the half-filled Hubbard model with a cut-off 
forbidding momenta close to the angles of the square shaped Fermi sur- 
face. By Renormalization Group methods we find a convergent expan- 
sion for the Schwinger function up to exponentially small temperatures. 
We prove that the system is not a Fermi liquid, but on the contrary it 
behaves like a Marginal Fermi liquid, a behaviour observed in the normal 
phase of high T c superconductors. 

1. Main results 



1.1 Motivations. The notion of Fermi liquids, introduced by Landau, refers to a wide class 
of interacting fermionic systems whose thermodynamic properties (like the specific heat 
or the resistivity) are qualitatively the same of a gas of non interacting fermions. While 
there is an enormous number of metals having Fermi liquid behaviour, in recent times 
new materials has been found whose properties are qualitatively different. In particular 
the high-temperature superconducting materials (so anisotropic to be considered essentially 
bidimensional) in their normal phase have a non Fermi liquid behaviour, in striking contrast 
with previously known superconductors, which are Fermi liquids above the critical temper- 
ature. While in Fermi liquids the wave function renormalization Z is Z = 1 + 0(A 2 ), where 
A is the strength of the interaction, in such metals it was found Z ~ 1 + 0(A 2 logT) for 
temperatures T above the critical temperature, see [VLSAR] (see also [VNS] for a review); 
metals behaving in this way were called Marginal Fermi liquids. Such results stimulated 
an intense theoretical research. It was found by a perturbative analysis, see for instance 
[AGD] or [Sh], that in a system of weakly interacting fermions in d = 2 Z is essentially 
temperature independent, at least for circular or "almost" circular Fermi surfaces. Despite 
doubts appeared about the reliability of results obtained by perturbative expansions [A], 
such results were indeed confirmed recently by rigorous Renormalization group methods. 
It was proved in [FMRT] and [DR] that indeed a weakly interacting Fermi system with a 
circular Fermi surface is a Fermi liquid, up to exponentially small temperatures. Such result 
was extended in [BGM] to all possible weakly interacting d = 2 fermionic systems with sym- 
metric, smooth and convex Fermi surfaces, up to exponentially small temperatures. These 
results cannot be obtained by dimensional power counting arguments as such arguments give 
a bound \Z — 1| < CA 2 | logT| from which one cannot distinguish Fermi or non Fermi liquid 
behaviour; for obtaining Z = 1 + (3(A 2 ) one has instead to use delicate volume improvements 



in the integrals expressing Z, based on the geometrical constraints to which the momenta 
close to the Fermi surface (assumed convex, regular and symmetric) are subjected. 

As Fermi liquid behaviour is found in systems with symmetric, smooth and convex Fermi 
surfaces, in order to find non Fermi liquid behaviour one has to relax some of such conditions. 
It was pointed out, see for instance [VR] and [ZYD], that the presence in the Fermi surface 
of flat regions in opposite sides could produce a non Fermi liquid behaviour; flat regions 
are indeed present in the Fermi surfaces of high T c superconductors [S] . The simplest model 
exhibiting a Fermi surface with flat pieces is the half-filled Hubbard model, describing a 
system of spinning d = 2 fermions with local interaction and dispersion relation given by 
e{k x , k y ) — cosk x + cos k v . The Fermi surface is the set of momenta such that e(k x , k y ) = 
and it is a square with corners (±7r, 0) and (0, ±tt). However this model has the complicating 
feature of vanishing Fermi velocity at the points (±7r, 0) and (0, ±7r) i.e. at the corners of 
the Fermi surface; this originates to the so called Van Hove singularities in the density 
of states. In order to investigate the possible non Fermi liquid behaviour of interacting 
fermions with a Fermi surface with flat pieces, independently from the presence of Van 
Hove singularities, one can introduce in the half filled Hubbard model a cut-off forbidding 
momenta near the corners of the Fermi surface. The half filled Hubbard model with cut-off 
(or the essentially equivalent, but slightly simpler, problem of fermions with the linearized 
dispersion relation e(k x , k y ) = \k x \ + \k y \ — it) has been extensively studied in literature, see 
for instance [M], [L], [ZYD] , [VR] , [FSW] , [DAD], [FSL]. The cut-off is somewhat artificially 
introduced but the idea is that the model, at least for same values of the parameters, belongs 
to the same university class of models with "almost" squared and smooth Fermi surface, like 
the anisotropic Hubbard models [Sh] , the Hubbard model with nearest and next to nearest 
neighbor interaction [M], or the half- filled Hubbard model close to half filling. 

Aim of this paper is to compute in a rigorous way the asymptotic behaviour of the 
Schwinger functions of the half filled Hubbard model with cut-off up to exponentially small 
temperatures. We will show that such a system is indeed a Marginal Fermi liquid, and our 
result furnishes indeed the first example rigorously established of such behaviour in d = 2. 

For our convenience, we will consider new variables k + = fcx + fc « and fc_ = fcx ~ fc " so that 
the dispersion relation of the half-filled Hubbard model is given by 

e(k + , fe_) = 2 cosfc + cos fc_ (1.1) 

and the Fermi surface is the set k + = ±-| or fc_ = ±f ■ 

1.2 The model. Given a square [0, L} 2 e R 2 , the inverse temperature (3 and the (large) 
integer M, we introduce in A = [0, L] 2 x [0, (3} a lattice Am, whose sites are given by the space- 
time points x = (xo, x+, x-) with (x+,X-) G Z 2 and xq = no/3/M, no = 0, 1, . . . , M— 1. We 
also consider the set T> of space-time momenta k = (fco, fc+, k-) = (ko, k), with k±, = 2 ™ ± , 
(n+,n_) € Z 2 , [-L/2] < n± < [L-l/2]; k = ^(n + |), n = 0, 1, . . . , M - 1. With each 

k e V we associate four Grassmanian variables ^ s , e, s e {+, — }; s is the spin. The lattice 
Am is introduced only for technical reasons so that the number of Grassmann variables is 
finite, and eventually the (essentially trivial) limit M — > oo is taken. Wc introduce also a 
linear functional P(dip) on the Grassmanian algebra generated by the variables V> k a , such 
that 

/ PW)K,sA,s* = L 2 l35 Sl ^ lM2 g{^) , (1.2) 
where g(k) is defined by 

m = *(k) (L3) 

— ik + 2 cos k + cos /c_ 

where x(k) is a cut-off function 

X (k) = H(a 2 sin 2 fc+)C Q - 1 (k) + H(a 2 sin 2 ft-JC^k) (1.4) 
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where 

C^ik) = H{^kl + 4cos 2 (fc+)cos 2 (fc_)) (1.5) 
and, if 7 > 1 and a > \/2 

The function C^" 1 (k) acts as an ultraviolet cut-off forcing the momenta k to be not too far 
from the Fermi surface, and ho not too large; the cut-off on ko is imposed only for technical 
convenience and it could be easily removed. The functions H{a^ sin 2 k±) forbids momenta 
near the corners of the Fermi surface i.e. the points (±7r/2, ±7r/2). The Grassmanian field 
V>x is defined by 

^ kex> 

The "Gaussian measure" P(dtp) has a simple representation in terms of the "Lebesgue 
Grassmanian measure" 

d^= n d$+ a dfo >a , (i.8) 

kei>,s=± 

defined as the linear functional on the Grassmanian algebra, such that, given a mono- 
mial Q(ip~ ,tp + ) in the variables tp^ , tp^ , its value is 0, except in the case Q(?p~ = 
riks^k spies' U P t° a permutation of the variable, in which case its value is 1. Finally 
Ilke-D s=± mcans a product over the k such that x(k) > 0. We define 

1 

P(#) =7V- 1 D^- exp[-— ]T X _1 (k)H/co + 2cosfc + cosfc_)^+^-J , (1.9) 

^ ker>,<T=± 

with TV is a renormalization constant and again ^ k means a sum over k such that x(k) > 0. 
The two point Schwinger function is defined by the following Grassman functional integral 

- y ) = J™ JP W )e-VM ' (L10) 
where, if we use / dx as a shorthand for -j^ X) x eA M ' 

V(V) = / ^V+, s V->+_>-_ s , (1.11) 
We call S(k) the Fourier transform of S(x — y). 

1.3 Main Theorem. Our main results are summarized by the following Theorem, which 
will be proved in the following sections. 

Theorem. Given a large enough, there exist two positive constants e and c such that, for 
all |A| < e an d T > cxpj-^IAI)- 1 }, for alike V such that ^ < ^ 7c 2 + 4 cos 2 k + cos 2 fc_ < 
and H(oq sin 2 fc_) = 1 then 

S(k) = (fco 2 +4cos 2 fc + cos 2 fc,)^-) A2 
v ; -iko + 2 cos k+ cos fc_ v w; v ; 

and /or keP such that ^ < ^//c 2 + 4 cos 2 k + cos 2 fc_ ) < || and i? (a 2 , sin 2 fc + ) = 1 i/ien 

S(k) = (fco 2 +4co S 2 fc + cos 2 fc-)^ + ) a2 
v ; -ifco + 2 cos k+ cos v v " v 7 
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where |A(k)| < c, where c > is a constant, and f](k±) = a(k±)X 2 + 0(X 3 ) is a critical 
index expressed by a convergent series with a(k±) > a not identically vanishing smooth 
function. 

1.4 Remarks. The above theorem describes the behaviour of the two point Schwingcr 
function up to exponentially small temperatures, i.e. T > exp{— (c| A | ) 1 } ; the constant c is 
essentially given by the second order terms of the perturbative expansion. A straightforward 
consequence of (1.12), (1.13) is that that the wave function renormalization is Z = 1 + 
0(A 2 log/3), which means that the half-filled Hubbard model with cut-off is a marginal 
Fermi liquid up to exponentially small temperatures. From (1.12), (1.13) we see that the 
behaviour of the Schwinger function close to the Fermi surface is anomalous and described by 
critical indices which are functions of the projection of the momentum on the Fermi surface. 
Critical indices which are momentum dependent were found for the same model also in 
[FSL] by heuristic bosonization methods. The presence of the critical indices makes the 
Schwinger function quite similar to the one for d = 1 interacting spinless fcrmionic systems, 
characterized by Luttinger liquid behaviour (see for instance [A]). However an important 
difference is that the critical exponent r\ in a Luttinger liquid is a number, while here is a 
function of the momenta. Another crucial difference is that in a Luttinger liquid S(k) ~ 
.g(k)|k| T ', with r\ = a\ 2 + 0(X 3 ) up to T = 0; hence a Luttinger liquid is a Marginal Fermi 
liquid for high enough temperatures but not all the marginal Fermi liquids are Luttinger 
liquids. 

The paper is organized in the following way. In §2 we implement Renormalization Group 
ideas by writing the Grassman integration in (1.10) as the product of many integrations at 
different scales. The integration of a single scale leads to new effective interactions, and the 
renormalization consists in subtracting from the kernels of the effective interaction, which 
are not dimcnsionally irrelevant, of the effective interaction their value computed at the 
Fermi surface. One obtains an expansion for the Schwinger functions as power series of a 
set of running couplings functions (depending from the momentum on the Fermi surface and 
the scale). In §3 we prove that this series is convergent if the running coupling functions are 
small enough; the convergence radius is finite and temperature independent, and this means 
that the theory is renormalizable. In the proof of convergence one uses the Gram-Hadamard 
inequality. In §4 we show that the running coupling functions obey to a recursive set of 
integral equation, called Beta function, and we show that the running coupling functions 
remain small up to exponentially small temperatures T > exp{— (c| Al)^ 1 }. Moreover we 
show that the wave function renormalization has an anomalous flow, with a non vanishing 
exponent (contrary to what happens for instance in the case of circular Fermi surfaces), 
and this essentially concludes the proof of the Theorem. It would be possible to use our 
beta function to detect (at least numerically) the main instabilities of the system at very 
low temperatures. At the moment, this kind of numerical analysis was done for this model 
only in [ZYD] in the parquet approximations, with no control on higher orders which are 
simply neglected. Finally in §5 we compare the Marginal Fermi liquid behaviour we find in 
this model with the Luttinger liquid behaviour, and we discuss briefly what happens in the 
Hubbard model with cut-off close to half filling. 

It is very likely that the half-filled Hubbard model with cut-off can work as a paradigm 
for a large class of systems, in which the Fermi surface is flat or almost flat but there are 
no Van Hove singularities. Marginal Fermi liquid behaviour can be surely found in the 
Hubbard model with cut-off and close to half-filling, up to temperatures above the inverse 
of the radius of curvature of the Fermi surface. Another model in which one could possibly 
find Marginal Fermi liquid behaviour is the anisotropic Hubbard model introduced in [Sh] 
with dispersion relation cosfci + tcos/c2, with t = 1 + e. Such model has a Fermi surface 
with no van Hove singularities and four "almost" flat and parallel pieces, and one can expect 
Z = 1+0(A 2 | log(|e|)| log/3) for (3 < 0(min[e _1 , exp{(c|A|)}]. Another interesting question is 
the possibility of Marginal Fermi liquid behaviour in the Hubbard model close to half-filling 
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(with no cut-off). At half- filling it is believed Z ~ f + 0(A 2 log 2 /3), so a different behaviour 
with respect to Marginal Fermi liquid behaviour. A renormalization group analysis for this 
problem was begun in [R], and it was proved the convergence of the series not containing 
subgraphs with two external lines for T > exp{ — (co|A|)~ }. 

2. Renormalization Group analysis 



2.1 The scale decomposition. As the spin index will play no role in the following analysis 
(on the contrary it is expected to have an important role at lower temperatures) we simply 
omit it. The cut-off function x(k) defined in (f .4) has a support in the k space which is given 
by four disconnected regions, each one containing only one flat side of the Fermi surface. 
It is natural then to write each Grassman variable as a sum of four independent Grassman 
variables, with momentum k having value in one of the four disconnected regions; each field 
will be labeled by a couple of indices, u = I, II and w = ±1, so that each field has spatial 
momenta with values in the region containing (up F ,0) if a = I or (0,u;p F ) if a = II. We 
write then Grassman integration as 

f p(#)f(v) = / n II p °A*i>)n E E ^.«) ( 2A ) 

J a=I,II u>=±l a=I, 1 1 w=±l 

where F is any monomial, u = ±1 and 

PiMUT _ $t /lt ,.,- =5 uu .&h k ,g( a 2 sin 2 fc_) C * Hfr), 

I,u,\ YJY IiU y i+upFtI Y IlU > lk >+u,'p FjI k p k v o J _ iko + 2ujsmk' + cosk- 

(2.2) 



/ 



/ 



C- 1 (k 07 k'_,k + ) 
iko + 2ui sin k'_ cos k+ 



r r tti 1 ■ 2 i \ C u 1 (fco J k'_, k + ) ,„,, 

= <5ki,k'5 w , U 'ir(aoSin K) (2.3) 



where 

Cj^fco, *+,*;_) =6{u>k' + + PF )H(^Jk% +4sin 2 fc^ cosfc_) (2.4) 

Cj^ifeb, #_,*;+) = e{ojk'_+p F )H{^kl+Asin 2 k'_cosk+) (2.5) 

and pV,cr is defined such that p F j = (f,0) and p F ,ii = (0, § ); moreover p F = | and 
k = k' + ijjp F , a ( fc' is the momentum measured from the Fermi surface). 

It is convenient, for clarity reasons, to start by studying the "free energy" of the model, 
defined as 

log J P{di>)e- V (2.6) 
where, calling with a slight abuse of notation Vv, w ,k'+u;pF a = Vv^.k', V is equal to 



r 4 

A E E /dki...dk^(^ £l (k^+^ CTi ))^+^ !k , 



^ (Jn .Uln .k' ^fT^ .QJq .kl 



CT 2 ,o; 2 ,k2 ^o-3,cj 3 ,k3 ^4,^4,^ 

ct?i,..,a?4 (7i ,..,<J4=/,// ^ i= 1 

(2.7) 

where / dk = ^ £ k and <5( k - k ') = L 2 (3S K w. 

We will evaluate the Grassman integral (2.6) by a multiscale analysis based on (Wilsonian) 
Renormalization Group ideas. The starting point is the following decomposition of the cut- 
off functions (2.4), (2.5) 

I : / : 

H(^fc 2 +4cos 2 Lsin 2 fc^ = ]T / fc (Vfc 2 +4cos 2 fc CT sin 2 fc;)= ^ fk{kvXX) 

A;— -co k— — oo 

(2.8) 
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with fk(t) = H{~f~ k t) — H(j~ k+1 t) is a smooth compact support function, with support 
7* : ~ 1 < \t\ < j k+1 ; moreover: 

&)k a = k- if a = I and k a = k + if a = II; k a is the projection of k in the direction 
parallel to the Fermi surface. 

b)fc^ = k' + if o = I and k' = k'_ if a = II; k' + ujpF,a is the projection of k in the 
direction normal to the Fermi surface. 

For each a, the function fk(ko,k!a, k a ) has a support in two regions of thickness 0{~f h ) 
around each fiat side of the Fermi surface, at a distance 0{^ h ) from it. We will assume 
L = oo for simplicity and it follows that there is a hp = 0(log/3) such that = for 
k < hp, while fk is not identically vanishing for k> hp. 

The integration of (2.6) will be done iteratively integrating out the fields with momenta 
closer and closer to the Fermi surface. We will prove by induction that it is possible to define 
a sequence of functions Zh(k' au} ) and a sequence of effective potentials V^ h ' such that 

J P / (#)^/(#)e- V = e~ L ^ J Pz h ,/W ( ^ ) )Pz h ,//W(^))e^ (h, (^ ( ^»). (2 .9) 
where Eh is a constant and y/Zh^-^ equal to 

(^Z h (k' IA )^%, ^ZhCk'^4-%, yJz h &iMw> \l Z ^'n,-^u%M>) ( 2 - 10 ) 
and Pz h ,a(dip(- h ^) is the fermionic integration with propagator 
, 1 H (a 2 . sin 2 k a )Cr},(ko, k' , k a ) 

) = 7 r v \ , — o * V *4 + "Pf) (2.11) 

Zj h\K a tj j) —iko + 2uj cos feo- sin fe^ 

with 

h 

C^JkoXX)= E h(k XX)- (2.12) 

/c— — oo 

The (^-function in (2.11) can be omitted by the definition of the variables k a . 

We define k„^ = (jj,ujpp,k-) if a = I and k' aul — (j|, fc + , wp^) if c = II; moreover 
k' a u = (|p 0, fc_) if (T = / and fc^ w = (-|, fc + , 0) if <r = 77; moreover we call k" u = (— -|, 0, fc_) 
if a = I and fc£ w = (-|, fc+,0) if a = II. 



If e = ± 



v^ ( ^) = E E E E 

n=l Ul,...,U2n ffl,...,J2» £!,...,£. 



^(k'l,-,^) (2.13) 



where 

^(kL.-kan-i) - Wt(K + WlPF, CTl ...k 2 „-l + Wan-lAr.^.J = ^(^...k^) (2.14) 



2.2 Tfte renormalization procedure. Let us show that (2.9) is true for h — 1, assuming 
that it is true for ft. We define an C operator acting linearly on the kernels of the effective 
potential (2.13): 

1) CW^ = if n > 2 

2) If n = 1 

£W 2 h (k>) = \[W 2 h (Kj + W$KJ] + k d ko W 2 h (K,J + smk' a dM(K,u,)} (2-15) 
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where dk means the discrete derivative and <9 CT = dk + is o = I and d a = dk_ is o — II. We 
will prove in §4 that [W^(~K,J + MiKJ] = 0. 

3)If n = 2 

£W£(k' l7 k' 2 X) = W^K^X^K^J (2-16) 

Calling d W£(k' u a ) = -ia^fc^J, d a W£(k' ut7 ) = 2wcosfc CT z/ l (A4 CT ) and 

lh(k' ai . lJl ,k' cr2 ^ 2 ,k' (J3 U)3 ) = W 4 (k' rTl ^ 1 ,k' rT2 0J2 ,k' t73 ^ 3 ) (2-17) 

we can write 

£V h = £ / ^(02^08^ sin ^ 

cr=IJI J M,M 

(2.18) 

We write the r.h.s. of (2.9) as 

J Pi,z h (d^ h) ) J P // ^J^^' l )) e - £v(h, (^ ( - h) )- KV< ' l) ^ ( - h) ) (2.19) 
with TZ = 1-C. 

2.3 Remark 1. The non trivial action of TZ on the kernel with n = 2 can be written as 

+ [^(^ 1;Wl ,k' 2 ,k^)-H^(fc; iia;i ,fc; 2;W2 ,k'3)] (2.20) 

+ [^4(^,^1 ^CT 2 ,0J 2 ' ~ ^4 (^CTi.O)!' ^CT 2 ,0) 2 ' ^3,0)3)] 

The first addend can be written as, if ai = I (say), in the limit L — > 00 

(fco,i ~ |) jf dtflfco,! + f(*o,i - ^,i> fc-,i; k 2' k s) + (2-21) 

fcV ;1 J dtd k , + W£(0, tk' + ^ k-,i; k' 2) k' 3 ) 

The factors fco,i — 7r//3 and fc^ x are 0(7^ ), for the compact support properties of the 
propagator associated to V^Lfk' ' w ^ n ^' — ^' wnue the derivatives are dimensionally 
0(7~ h+1 ); hence the effect of TZ is to produce a factor 7^ -k^ 1 < 1. Similar considerations 
can be done for the other addenda and for the action of TZ on the n = 1 terms. 

Remark 2. From (2.16) we see that the effect of the C operation is to replace in W^ik) 
the momentum k with its projection on the closest flat side of the Fermi surface. Hence 
the fact that the propagator is singular over an extended region (the Fermi surface) and not 
simply in a point has the effect that the renormalization point cannot be fixed but it must 
be left moving on the Fermi surface. 

2.4 The anomalous integration. In order to integrate the field ipW we can write 

J P JiZfc ( d ^ fc >) J P/i,z k (# ( ^ ) )e- £v(k, (^ (Sh, »- RvW ^ (Sk, » = (2.22) 
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I P IiZh _Ad^) I P // ^ h _ i (^^' l )) e - £1>h (^<^ ) )-Kv('->(^<^') 

where P a ^z h _ 1 (dip(- h ' 1 ) is the fermionic integration with propagator 

1 Hjals^k^C-^kvXX) (2 23) 
Zh-i(k') — ik Q _|_ 2u>cosk a sinfc^ 

and 

Z h _!(k') -^(fc^)[l + ff(agsin 2 fc CT )C ft - 1 (fco,fc;,fc CT )^ i (fc^)] (2.24) 

Moreover 

CV h = CV h - I dtiz h {k^ a )[2u C osk a sinfc; - ik ]^ h J^ h J. (2.25) 

We rescale the fields by rewriting the r.h.s. of (2.22) as 

J Pi^Ad^ h) ) J P /7 ^ h _ iW (^' l ))e-^ h (V^ 77 ^ < - h) )- KV<h) (V^ 77 ^ ( -' l) ) (2.26) 
where 

CV h = fdkJ2 [ShAK)^™sKsmgJ]^J- M ,^+ £ [dK-dK 

J <y=I,II a 1 ,..,a i =I,II'' 

^(^^^k^^ (2.27) 

i 

and 

Sh(K,a) = - *h(K,j) (2.28) 

Z h-l(k') 



4 



We will call <5^ and A/j running coupling functions] the above procedure allow to write a 
recursive equation for them, see §5. 
Then we write 

J P I , Zh _ 1 (d^- h - 1) ) J Pii^Ad^ h - iy ) J Pi,z h _Ad^ h) ) J Pn^Ad^ h) ) 

and the propagator of P rT ^z h _ 1 {dtp) is 

gtA*)=H(al^k a ) L AfoX.M (2.30) 

z 'i-il fc w,crJ — «fco + 2uj cos smfc CT 

and 

with -ff(ao sin 2 fc CT )/^(fc , fc^, fc CT ) having the same support that H(clq sin 2 k a )fh(ko, fc^i 
We integrate then the field i/'' 1 and we get 



e -L J /3E h 



1 / P/,z,_ 1 (# ( ^- 1) ) / Pn, ZH (# ( ^ i y vM,( ^"" ,) (2.32) 
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and the procedure can be iterated. 
We will see in the following section that if the running coupling functions are small 

sup sup \6 k (K,J\ < 2|A| sup sup ^''^ < e 2|A| 

sup sup |A fe (^ lia)1 ,^ 2 , W2 ,fc; 3ia) 3)|<2|A| (2.33) 

then the effective potential is given by a convergent series. In §4 we will show that up to 
exponentially small temperatures this is indeed true. 

3. Analyticity of the effective potential 



3.1 Coordinate representation. It is convenient to perform bounds to introduce the vari- 
ables V'xoi ct- We define the fields = e ieul ^ F '" x ^%_ IJ1 or more explicitly 

= e*""""^,! r XtU ,n = e tEUJPFX2 t>Un (3-1) 
and the propagators of such fields is 

t M - y) = / dk' 1 e-^lx-y) H ^ Sin2 VMkoXX) (3 2) 

J Zh-\(k' u a ) — ik + lo2 sin k a cos k a 

It is easy to prove, by integration by parts, that for any integer N, for L — > oo 

C ^(l+Il +ll + ) 

IWt^U* - y)i < l + m ^^L-^-^r <») 



(3.4) 



' ' ' " l + h h \d(x a -y )\ + \x + -y + \+^\x_-y_\r 

where d(x ) = § sin 

Proof. The above formula can be derived by integration by parts; note that, if for instance 
a = I 

<9fc_ — ; ^ — : — n ; — = 7 — ; 7. — : — n ; — r^"2w sin k' + sin fc_ (3.5) 

— iko + 2u sin k + cos k- {— ik n + 2uismk + cos/c_)^ T 

which is 0(j~ h ); in the same way the n-th derivative with respect to fc_ is still 0(j~ h ). On 
the other hand d%°d%+ is bounded by ^-h-n Q h-n + h.^ nna Uy the integration gives a volume 
factor 7 2h . 

We define 

J<( Xl ,...,x 2 „) = -i-^ ^ e- < S"i e ^T^(ki,...l4 l _ 1 )«yE e< (k{+a; i ft,, <ri )) 

1 J ki.-,k 2 „ 

(3.6) 

Hence (2.13) can be written as 



vW(^ fc ) = E E E E /^...dx 2 



n = l Wl,...,U2n fi,...,<T2n £l,---,£2n ' 



2n 



W^^Xi, ...,x 2 „) 

(3.7) 
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We now discuss the action of the operator £ and 1Z = 1 — C on the effective potential in 
the re-space representation. Noting that from (3.6), if £1 = £ 2 = —£3 = —£4 = + 

iy 4 /l (x 1 ,x 2 ,X3,x 4 ) = e" 4 ^ 1 ^^ 2 ^-^,^-^^ (3.8) 
we can write the action of 1Z (2.16) as 

4 4 

J i=l i=l 

/: 



= /rr^ xi n^' ei e ix 4 (wiPF, CT1 +w2PF,<T 2 -w3PF, CT3 -w 4 p Fi<T4 ) 
i=i i=i 

[W4 (xi - x 4 ,x 2 - x 4 ,x 2 - x 4 ) - S(x ,i - x a)S{x .2 - x , 4)6(x , 3 - x 0A )5(x ail - x^ lA ) 

)5(^t 3i 3-^3,4) / ^o,i^o,2rfto,3<iii >ffl rft 2iCT2 (ii3 iCT3 W 4 ( ' l) (ti,t2,t3) (3.9) 

where / efoc = X) xeA and *» = ^0,1,^,^,1) where i J;i = = t-,i and = 

On the other hand we can equivalently write the 1Z operation as acting on the fields, and 
such two representations of the 1Z operation will be used in the following. It holds that, by 
simply integrating the deltas in (3.9) 

4 4 

J 1=1 J 1=1 

\D + - i> + ib~ i>~ +H>± D + - tb~ i)~ + 

L Xi,X4 I „ 1 , <7l,Wl ^X 2 , <T2,^2 r X 3 ,(T3,W3 Vx 4 ,<T4,UJ 4 1 rx 4 ,CTl ,Ui X 2 ,X 4> „ 2 ,CT 2 ,W2 ^X 3 ,<J 3 ,W 3 V x 4 ,<7 4 ,W 4 

(3.10) 

where x 4 ^ CTi = (#0,4, sc+,4, if cr, = / and x 4j(Ti = (x ,4, ^-,4) if o» = ^> moreover 

^XdiU.jj^dUi = V'xjjO-jjWj — V'I 4i0 . i ,0-4,0)4 (3-11) 

This means that the action of the renormalization operator 7?. can be seen as the replacement 
of a xp e( -- h ~> field with a D^^ a field and some of the other ^'(-'O fields are "translated" 

in the localization point. The field D^r^^ ff Wj is antiperiodic in the time components of 
x i; and x 4 . We can write D £ (<' 1 ) as sum of two terms (if at = I for instance): 

^Xi,X4 j „ i ,cr i ,a>i [V'xi.j Vxo^jX+^x-,,,/] "I" [V'xo^jX+.^x-,,,/ V'xo^jX+^x-,,,/] (^-12) 

and the second addend can be written as, for L — > 00 

V'x 0i 4,x + ,i,x_,i,J ~~ V'xo, 4 ,x + , 4 ,x_ j ,,7 = ( x +,i ~ ^+,4) / ^X + '^/, w ,x ,4,X + ,i-t(x + ,i-X + ,4),X_, i 

J 

(3.13) 

and x + ,i — t(x+ t i — x +A ) = x + ,i A (t) is called interpolated point. 

This means that it is dimensionally equivalent to the product of the "zero" (x + .i—x +A ) and 
the derivative of the field, so that the bound of its contraction with another field variable on 
a scale bl < h will produce a "gain" j~( h ~ h ), see (3. 3), (3. 4), with respect to the contraction 
of ip l 0? >a ■ Similar considerations can be repeated for the first addend of (3.3); some care 
has to be done as (3 is finite, and we refer §3.5 of [BM], If there are two external lines. 



TZ 



J dxidx 2 W 2 ,l (xi,X2)^+ ) 



<Ti,Wi Vx 2 ,<T 2 ,w 2 



x 2 ,Xi >CT2 ,cr 2 ,o) 2 



(3.14) 
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where 

T x 2 ,i! = ^,^,0,-^1,^,0,,^ - 0*0,2 - Zo,i)<9o#i „ 2 , W , CT2 - {x aa -x lta )d a i>^ a2>w>(72 

(3.15) 

and ^ = 9^ if a = I and 9 CT = d x _ if <r = II. In this case the "gain" produced by the 1Z 



operation is 7 



-2(h-h') 



We can write the local part of the effective potential (2.18) in the following way 
CV h = J ^^[^((^-^Je^ 



(3.16) 



^ / dX4dXv 1 ,idx a2 , 2 dx c73 ^\h-, ull ,..oj 4 ((xi - £4)0-1, (^2 - Xi) ai , (X3 - £4)03) 

CTl,..,CT 4 =i",// 

.^(wipjf.oj+^Pf.oj -U3W , ,»3-W4PF,<, 4 )I+ ,7,+ ,7,- ,7- 

r X4 i „ 1 ,cri ,0,1 r x 4j0 - 2 ,02,0,2 t / x 4j0 - 3 ,0-3,0,3 rx4,a 4 ,U4 

where 5h(x) is the Fourier transform of <5/ l (fc cr )2 cos k a with respect to k a and (2^ — Xj) a — 
x-.i — X-j if a = I and (xi — £j) CT = a; +i , — x + j if a = II; moreover x 0i = X- if a = I and 
x ai = x + if a = II. 

3.2 Tree expansion. By using iteratively the "single scale expansion" we can write the 
effective potential V^^-^), for h < 0, in terms of a tree expansion. For a tutorial 
introduction to the tree formalism we will refer to the review [GM]. 




We need some definitions and notations. 

1) Let us consider the family of trees which can be constructed by joining a point r, the 
root, with an ordered set of n > 1 points, the endpoints of the unlabeled tree (see Fig. 1), 
so that r is not a branching point, n will be called the order of the unlabeled tree and the 
branching points will be called the non trivial vertices. The unlabeled trees arc partially 
ordered from the root to the endpoints in the natural way; we shall use the symbol < to 
denote the partial order. 

Two unlabeled trees are identified if they can be superposed by a suitable continuous 
deformation, so that the endpoints with the same index coincide. It is then easy to see that 
the number of unlabeled trees with n end-points is bounded by 4". 

We shall consider also the labeled trees (to be called simply trees in the following); they 
arc defined by associating some labels with the unlabeled trees, as explained in the following 
items. 

2) We associate a label h < — 1 with the root and we denote Tu,n the corresponding set of 
labeled trees with n endpoints. Moreover, we introduce a family of vertical lines, labeled 
by an integer taking values in [h, 1], and we represent any tree r € Th,n so that, if v is an 
endpoint or a non trivial vertex, it is contained in a vertical line with index h v > h, to be 
called the scale of v, while the root is on the line with index h. There is the constraint that, 
if v is an endpoint, h v > h + 1. 
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The tree will intersect in general the vertical lines in set of points different from the root, 
the endpoints and the non trivial vertices; these points will be called trivial vertices. The 
set of the vertices of r will be the union of the endpoints, the trivial vertices and the non 
trivial vertices. Note that, if v\ and v 2 are two vertices and v\ < v 2 , then h Vl < h V2 . We 
will call s v the number of subtrees coming out from v. 

Moreover, there is only one vertex immediately following the root, which will be denoted 
v and can not be an endpoint; its scale is ft + 1. 

3) To each end-point of scale +1 we associate V (1.11). With each endpoint v of scale 
h v < we associate one of the two terms appearing in (3.16), with coupling Xh v -i or 5h v -i- 
Moreover, we impose the constraint that, if v is an endpoint and h v <0,h v = h v > + 1, if v' 
is the non trivial vertex immediately preceding v. 

4) We introduce a field label f to distinguish the field variables appearing in the terms 
associated with the endpoints as in item 3); the set of field labels associated with the 
endpoint v will be called I v . Analogously, if v is not an endpoint, we shall call the set 
of field labels associated with the endpoints following the vertex v; x(/), e(f) and oj(f) will 
denote the space-time point, the e index and the u> index, respectively, of the field variable 
with label /. 

If h v < 0, one of the field variables belonging to I v carries also a derivative d a if the 
corresponding local term is of type S, see (3.16). Hence we can associate with each field 
label / an integer m(f) € {0, 1}, denoting the order of the derivative. 

If ft < — 1, the effective potential can be written in the following way: 

oo 

V^(^) + L(3E h+1 = J2 E VW(r,^ h) ), (3.17) 

where, if v is the first vertex of t and n, .., t s (s — s Vo ) are the subtrees of t with root v , 
V^ h \T,ip^- h ^) is defined inductively by the relation 

V("»(r,^f)) = 

(-1V+ 1 ( 3 - 18 ) 
and V {h+1 \n,^ h+i y) 

a) is equal to TZV^ h+1 \ri, ip(- h+1 ' 1 ) if the subtree Tj is not trivial, with 1Z defined as acting 
on kernels according to (3.9) and its analogous for n = 1; 

b) if n is trivial and h < -1, it is equal to CV {h+1 \^ h+ ^) (3.16) or, if ft = -1, to V. 

££ +1 denotes the truncated expectation with respect to the measure Yl a=I u Pa,z h (dip^ h+1 ^ ) , 
that is 



. (3.19) 

A,=0 



We write (3.18) in a more explicit way. If ft = — 1, the r.h.s. of (3.18) can be written in 
the following way. Given r G T_i^ n , there are n endpoints of scale 1 and only another one 
vertex, vq, of scale 0; let us call vi, . . . , v n the endpoints. We choose, in any set I Vi , a subset 
Q Vi and we define P Vo = DiQ Vi . We have 

y(-i)( T; ^(<-i)) = V^(t, P vo ) , (3.20) 

Pvc 

V^\t,P V0 ) = J d^J^iP^K^M , (3.21) 
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n 

K r,pJ^o ) = -jSo [V' (0) ( P V1 \Qv! ),•••, ^ (0) (P Vn \Qv n )] J] K V? ( X "* ) . ( 3 - 22 ) 

i=l 

where we use the definitions 

^)(P„) = J] ^^igf (/) , h<-l, (3.23) 

fePv 

&°hp v )= n^°(/) (/) ' ( 3 - 24 ) 



x-vi = x (3.25) 

It is not hard to sec that, by iterating the previous procedure, one gets for 1/( ,1 )(t,V> ( -' i) ), 
for any r G T^n, the representation described below. We associate with any vertex v of the 
tree a subset P v of 7„ , the external fields of u. These subsets must satisfy various constraints. 
First of all, if v is not an endpoint and v\ , . . . , v Sv are the vertices immediately following 
it, then P v C UiP„ 4 ; if u is an endpoint, P v — I v . We shall denote Q Vi the intersection 
of P v and P Vi ; this definition implies that P v = UiQ Vi . The subsets P Vi \Q Vi , whose union 
T v will be made, by definition, of the internal fields of v, have to be non empty, if s v > 1. 
Moreover, we associate with any f E T v & scale label h(f) = Given r € 7^ >T j, there are 
many possible choices of the subsets P v , v e r, compatible with all the constraints; we shall 
denote P T the family of all these choices and P the elements of V T . 
Then we can write 

V( fc )(T,V ( - fc) )= 51 V ( ' l) (r,P) • (3.26) 

yW(r, P) can be represented as 

VW(t,P) = | dx.^^' l )(P„ )^ p +1) (x„ ) , (3.27) 

with K^ h p ^(x„ ) defined inductively (recall that /i„ = h+ 1) by the equation, valid for any 
w G r which is not an endpoint, 

iffp } (x„) = J- JJ[^+i)(x„J] ^[^>(P t , 1 \Qt, 1 ),...,^ (fc - ) (i , „.„\Q w .J] , (3-28) 

where ^ hv \P v ) is defined as in (3.23), with (h v ) in place of (< h v ), if h v < —1, while, if 
h v = 0, it is defined as in (3.24). 
Moreover, if v is an endpoint and h v — 0, K^\yi v ) is given by (3.25), otherwise, see (3.16) 

K (h v ), x ^ = Uh v -i(xi,X2,x 3 ,x 4 ) if v is of type A, , g 2Q . 

" " \ rfft„-i(a f i, 2^2) if w is of type z, 

where 

e af4( E io; 1 ft ? ,„ 1 + ea ^ F .. a -a,,fe„3- e ^4?F..J A ^_ 1)(i;i) __ )(i)4( ( £l _ £ 4 ) CTl , (£ 2 - x 4 )„ 2 , (£3 ~ X*)a 3 ) 

If v is not an endpoint, 

Ki h ^\^ Vi )=nK^l^ Vi ) , (3.30) 

where Ti is the subtree of r starting from u and passing through i>, (hence with root the 
vertex immediately preceding v), PW and is the restrictions to of P. The action of 1Z is 
defined using the representation (3.9) of the 1Z operation. 

(3.26) is not the final form of our expansion, since we further decompose V^ h '(r, P), by 
using the following representation of the truncated expectation in the r.h.s. of (3.28). Let 
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us put s = s v , Pi = P Vi \Qvi', moreover we order in an arbitrary way the sets = 
{/ G Pi,a(f) = ±}, we call their elements and we define = Uy gP -x(/), yM = 
U /£P +x(/), XlJ = x(/r.), yij = x(/+.). Note that £? =1 I^TI = £i=i l^ + l = ™> otherwise 
the truncated expectation vanishes. A couple I = (fij,f^j,) = (/;~>/; + ) wm be called a 
line joining the fields with labels f^,f#j, and sector indices lo^ = w(fj^), = and 
connecting the points x; = Xjj and y; = y,'j', the endpoints of /. Moreover, we shall put 
to; = m(fj~) + to(/ ; + ) and, if o> ; ~ = = = A similar definition is repeated for 

a. Then, it is well known (see [Le], [BM], [GM] for example) that, up to a sign, if s > 1, 

swvi). -,j> w {p,)) = e n i- a - ^ - 6 °T,°t I dPT{t) dct G " ,T(t) 

(3.31) 

where T is a set of lines forming an anchored tree graph between the clusters of points 
X W u yW ; that is T is a set of lines, which becomes a tree graph if one identifies all the 
points in the same cluster. Moreover t = {ti^ G [0, 1], 1 < < s}, dPx{t) is a probability 
measure with support on a set of t such that t^i> = Ui ■ Uj' for some family of vectors e R s 
of unit norm. Finally G h,T (t) is a (n — s + 1) x (n — s + 1) matrix, whose elements are given 

by G^j, = tiyd^i^d™^^ with (/y./.tj-,) not belonging to 

T. 

In the following we shall use (3.31) even for s = 1, when T is empty, by interpreting the 
r.h.s. as equal to 1, if |Pi| = 0, otherwise as equal to detG' 1 = E^{^ h \P\)). 

If we apply the expansion (3.31) in each non trivial vertex of t, we get an expression of 
the form 

VW(r,P) = W ^^^(PJflJ^x.) = Y, V W (t,P,T) , (3.32) 

where T is a special family of graphs on the set of points x„ , obtained by putting together 
an anchored tree graph T v for each non trivial vertex v. Note that any graph T € T becomes 
a tree graph on x„ , if one identifies all the points in the sets x v , for any vertex v which is 
also an endpoint. 

We are writing the TZ operation as acting on the kernels, according to (3.9) and its analo- 
gous for n = 1 . Such representation for the TZ operation is however not suitable to " gain" the 
convergence factor j~( h ~ h ), or -y^ 2 ( h ^ h ) ; for which is much more convenient representation 
of TZ in (3.10), (3.14). However if we write simply all the TZ operations as in (3.10), (3.14) one 
gets possibly factors (xj — Xj) a ™ with a n — 0(n), which when integrated give 0{n\) terms. 
One has to proceed in a more subtle way starting from the vertices of r closest to the root 
from which the TZ operation is non trivial, and writing TZ as in (3. 10), (3. 14) leaving all the 
other TZ operation as in (3.9). One distributes the "zero" along a path connecting a family 
of end points, and from (3.9) (xj — Xj)TZW^ = (xj — Xj)W^, if Xj,Xj are two coordinates 
of W% and TZW^ is the term in square brakets in the l.h.s. of (3.9); an analogous property 
holds for TZW%. There are same technical complications in implementing this idea, which 
are discussed in [BM] (see also [BoM]), §3.2, §3.3 for a different model, but the adapting of 
such argument to the present case is straightforward. We obtain, in the L — > oo limit 



y^(r,p) = E E / d Xvo ^ lPvol w T , P , T , a (x vo ) ■ 

\ 11 ^a(/)^x a (/), W (/)/ ' 



feP v , 
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where 

APv\/2 



W, 



■,P,T,a(x„ ) = [ [] { Z hJ Z K- 



v not e.p. 



•[n4:K)( x -y^( x <)]{ II f dPrM • (3.34) 

i=l v not e.p. 11 ' 

Iet v 

where: 

1) P is the set of {P v }; 

2) T is the set of the tree graphs on x„ , obtained by putting together an anchored tree 
graph T v for each non trivial vertex v; 

3) At is a set of indices which allows to distinguish the different terms produced by the 
non trivial TZ operations and the iterative decomposition of the zeros; u*,.. are the 
endpoints of r, ff and / ; + are the labels of the two fields forming the line /, "e.p." is an 
abbreviation of "endpoint" . 

4) G^ l v ' Tv (t v ) is obtained from the matrix G h "< Tv (t„), associated with the vertex v and T v , 
by substituting G^Jj, = U^dT^-A U f + J '\g[ h - ) - + , ,+ (xy - y^) with 

5) 9? , g = 0, 1, 2, are discrete derivatives or operators dimensionally equivalent to derivatives, 
due to the presence of the lattice and the fact that (3 is finite, see [BM] §3. Morever 8® denotes 
the identity and j = 0,+,—. According to (3.13), (3.15) if <r(/) = / then in d?,~ one has 
j(f) = 0, + and if a(f) = II then j(f) = 0, -. 

6) d (x ( - yi) = | sin ^(x ,i - yo,i) and <i(x ( - y ; ) = (x it i - y it i), i = ± are the "zeros" 
produced by the TZ operation, see (3. 13), (3. 15). Finally by construction b a (l) < 2. 

7) The factors z ^~ 1 arc functions of the coordinates, and such dependence is not explicitly 
written. 

Of course the coefficients b a and q a are not independent, and, by the definition of TZ (see 
the discussion after (3.13)) it holds for any a G At, the following inequality 

[ Y[ 7 M/)M/)][I] 7 -M0M0] < Yl (3-36) 

/e/„ ZeT ii not e.p. 

where h a (f) — h Vo — 1 if / G P„ , otherwise it is the scale of the vertex where the field with 
label / is contracted; h a (l) — h v , if I G T v and 

fl if|P„|=4, 
= ^ 2 if |P„|=2; (3.37) 
I otherwise. 

It holds 

\detG h a ^(t v )\<C^Z^-^- 2 ^ ■ 

. 7 ^(E*=il p ««l-l^l- 2 ('«- 1 )) 7 fc «E;:i[«-( J '«*\««*)+ TO ( J '«*\«««)] . (3.38) 

. 7 -fc«E, 6 T [««"(A + )+9«.wr)+ m (/i + )+ ro wr)] . 

This follows from the well known Gram-Hadamard inequality, see also [Le] , [BM] , [GM] , 
stating that, if M is a square matrix with elements My of the form M y - =< Ai,Bj >, 
where Aj, _Bj are vectors in a Hilbert space with scalar product < •, • >, then 

|detM|< JjHAill-HBill. (3.39) 
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where || • || is the norm induced by the scalar product. 
In our case it can be shown that 



r h v ,T v _. ,n q <*( f i^n q <x( f &xrnVr)ftrn(f+)(h v ) (v \ 



= < 



(3.40) 



where e R s , i = 1, ...,s, are the vectors such that j/ = Ui • Uj/, and , 

B^ hv \ , are such that (in the case g = m = for simplicity): 

x (/,Tj/),"i,o-i 



O.w^o-i ' x(/+ ,),ui,ai 



= I dyA^/V B {K \ f+ 

' x(/. .)-y,o)(,cri y-x(/T ,),Wi,<7i 



(3.41) 



For instance A and _B can be chosen as: 
A?*] =-iJ dk' 



-ik'y 



H(al sin 2 k ai ) f h (k , ,'k ai )- 



2 + (2cosfc CT! sinfc'J 2 2 



B&lv, = J dk ' e ^'*\J H ( a l sin 2 fc CTi )/ft(fc ,^ , ifc + 2w cos sin fc' ; 



and from (3.39) we easily get (3.38). 
By using (3.34) and (3.38) we get, assuming (2.33) 

,ct( x u )l — C Jr.p,r,T,a j J 

j; not e.p. 

. cEii i^j-i p "I- 2 ( s "- i ) 7 t(e;: 1 1^1-1^1-2(^-1)) 



(3.43) 



• 7 



where 



/n 
rfx„ 1 [ JJ j (xi , yi ) AT^ (x„ . ) 
i=i 

■ { n U n ^ WHS^.^ir-.vf <»> - «»]} 



■unot e.p. ZeT„ 

In [BM],[BoM] it is proved that 



d(x t , ) = dx]Jdr; 



(3.44) 
(3.45) 



where r; = xj(t;) — y;(sj) and Xj(fj), yj(s;) are interpolated points, see (3. 13), (3. 15), and x 
is an arbitrary point of x„ . By using (3.3), (3.4) we bound dimensionally each propagator, 
each derivative and each zero and we find 



j TPTa <r TT [— c 2 ^-V^ Ei ^ MO+Mo - 

' ' ' LL IsJ 



d not e.p. 
JieT v 



. 7 -^( S ^i) 7 ^E, er J^(/ ! + )+^(/r)+™(/ i + )+™(/D] 



(3.46) 



We find then 



J dx„ |W T)P) T,a(x„ )| < 
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C n L 2 p\X\ n l- hD(P ^ II f^ T cS:: i l^l-l^l 7 -[-2+^+^)]| ( 3.4 7) 

v not e.p. 

where D(P V0 ) = — 2 + m. The sum over t, P, T, a is standard and we refer to [BM], §3.15; 
at the end the following theorem is proved. 

Theorem. Let h > hp > 0. If (2.33) holds then there exists a constant Co such that 

E E E E f dx V0 \w T , P ,T, a (*v )\ < 

TeTh -" ip va ^=2 m TeTaeAT (3-48) 
<L 2 fa-Woo) (co\) n , 

where 

D(P V0 ) = -2 + m. (3.49) 



4. The flow of running coupling functions 



.1 Lemma It holds that J2 V =± Wj Hv%, ±f ,*-) = Er,=± W i ! fajj ,k+, ±f ) = 



Proof - We can compute W^k) also by a "single scale" integration; in fact W^(k) is the 



kernel of the term ipu^-^^k on V '' l defined by 



e- v '' h = J P{d^)e- v ^° ]+ ^ <h) ) (4.1) 
where P(cfo/>'' 1 ' ') is the fermionic integration with propagator 

^ ( °' +' ) — ifc + 2 cos fc + cos fc_ ^ ^ 

where 



Xfc,o(k) = [i?(ao sin2fc +) + ^( a o sin2fc -)]^ 1 (\/^o +4cos 2 fc+cos 2 fc_) (4.3) 

with C^ 1 = Yl°k=hfk- We can write W^ijj, f, as sum over Fcynman graphs (see for 
instance [GM]) and each Feynman diagram can be written as 

7^E---7^E^k 1 )... 5 [ ^^ (4.4) 

ki k / 1 t — 1 i — 1 

where n + m is an odd number, a 1 , a\ = 0, 1, — 1, a 3 + J2i a l 1& an °dd integer and k = 
(fifi^)- In order to write (4.4) we consider a spanning tree T formed by propagators 
connecting all the vertices of the graphs. We will call the propagators not belonging to 
T loop lines and we write the momenta of the propagators of T as a linear combination 
of the momenta of the loop propagators and of the external momentum. We perform the 
shift k +: i — ► k' +i + §, and the summation domain is not changed by periodicity. The loop 
propagators become 

jM]( k ') = Xh,o( k ') (45) 

— ik + 2 sin k' + cos k- 
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with 

Xh.oflO = [ H ( a o cos2 k '+) + H ( a l sin2 k-)]C^(yJk^ + 4 sin 2 jfe^ cos 2 (4.6) 

Of course g^.o]^) is odd in the exchange /c ,fc+,fc- — > — fc ,— fc + ,fc_. On the other hand 
the momenta of the propagators belonging to T becomes 

m m ra ra m 

^^k J+ ^k = (^^fc , + ^^E^ fc ^ + (E^+^)f'E^ fc ^+ ff ^-) ( 4 - 7 ) 

i— 1 i— 1 i— 1 i— 1 i— 1 

with (X^i °i + <T ' 7 ) an odd integer; hence the propagators belonging to T have the form 

XhA k ') , 481 

-iEi + CT f ] + 2 K- 1 )' C sin(E, ^^, + )Hcos(Ei + *M ' 

and 

X M (k') - [H(al cos 2 (^ ct^^J + if (a 2 sin 2 (]T ai k- ti + fc_))] 

[C^{ l{J2 vM 2 + 4sin 2 (^ <7^ + ) cos 2 (^ <7<A:_,j + <rfc_)) (4.9) 

V i i i 

Hence by performing the change of variables k ,k+,k- — > — fc , — wc find 

^,|,*-) = -W J h (-^|,fc_) (4.10) 



4.2 Finite temperature flow. The multiscale analysis defined above has the effect that the 
running coupling functions 5 h {k' a u ), Z h {k' a (jJ ) and Ah(fc^ 1)a)1 , k' U2 ^ , fc^ 3)a)3 ) verify a recursive 
relation of the form 

/^A (^(71,0)1 ' ^<T 2 ,W2' ^3,W 3 ) 

It is quite easy to prove that, at temperature not too low, indeed (2.33) hold. The proof is 
done by induction assuming that (2.33) holds for h and proving that it holds also for h— 1, if 
h— 1 > /i/3 and /3 < expc|A| _1 , where c is a suitable constant. In fact iterating for instance 
the last of (4.11) we find 

o 



fc=/l+l 



and from (2.33) and (3.48) we find, for |A| < 



2c o 

I sup PiiK^X^X^ <2c 2 A A 2 

{fe'},{(T} 



if 4 > is a bound for the norm of the second order contribution to A^. Hence 

sup | A ft _ i (k ai iUl , k a2 )W2 , fc CT3 )a)3 , ^ CT4 )W4 ) | < [|A| + |^|2c2 A 2 ] . (4.13) 
{ fe },{ CT } 
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Then sup |Afc_i| < 2|A| if (3 < e 2c 2 |A| , as \h\ < \hp\. The same argument can be repeated for 
S h and and (2.33) holds. 

4.3 Flow of the wave function renormalization. To complete the proof of the main Theorem 
one has to check that indeed the critical index 7/(fc+) or rj(k-) are non identically vanishing, 
and this is equivalent to show that there exists a non vanishing function a(k' a w ) > such 
that 

e-^f^^h < z h {k' auJ ) < e -MK. u )^h ( 4 U ) 

From the fact that /^(fc^J = £~ =2 (K,J with |0? (n) (*£,JI < c^|A|™, as a conse- 
quence of (3.49) and (2.33), it is sufficient to find an upper and lower bound for (3^ 2 \ From 
an explicit computation one finds 

2a ) cosU3f i) =24 E ^"[/^i^dksffl^^i) ( 4 -!5) 



5a 2 ,^ 2 ( k 2)5^ 3 ( k 3)<5( k + k 3 - ki - k 2 )X h (k CT , k 3j(T3 , ki )(7l )A/ l (ki )(7l , k 2 , CT2 , k CT )] | k= ^ ^ . 

where = J2 k=h/} 9v,w As the dependence from the momenta of Xh is quite complex, 
it is convenient to replace in the above integral Xh with A; if the integral so obtained is 
nonvanishing, the correction will be surely smaller for T > e~( c l A " for a suitable c, as 
Xh = A + 0(A 2 1 log from (4.13). We can choose a = I for definitencss (the analysis 
for a = II is identical), and we can distinguish two kind of contributions in the sum over 
°~\ i °2 1 o~3 ! one in which all the propagators are <?/ , and the other such that there is at 
least a propagator g u . The estimate of this second contribution is 0(X 2/ y h ), as it can be 
immediately checked by dimensional considerations and applying the derivative in (4.15) 
over the gu propagators (one can always do that). We can further simplify the expression 
we have to compute noting that 

with 

'^ (X - 7)1 " ^ \ + b ^-yo\ + ^-y + \ + \ X .-y.\r ^ 

i.e. similar to (3.3) with an extra a$ 2 . We can replace in (4.15) the propagators g\ h l with the 
first addend in the r.h.s. of (4.16); if such term will be given by a nonvanishing constant, 
the correction will be surely smaller at least for ao large enough. Hence the dominant 
contribution to (4.15) is given by 

J dk^ A dk^^H(alsin 2 k^ A )H(alsin 2 k^^)H(alsm 2 (-k^ A + k^^ + k^))A (4.18) 



with 



A= ^2 J dk ,idk ,3 J dk +t idk +t3 - 



/^(fco,i,*:+,i) 



-ikn i + u>i2k + i 

W1,W2,W3 " 

/^(fc , 3 ,fc + , 3 ) fH-koA + fcp, 3 + fcp, -fc + ,i + fc+,3 + fc+) , _ _ , , 

-zfc 0i3 + ^ 3 2fc+, 3 k+ -*(-fcd,i + fco,3 + fc ) + 2co 2 (-fc +!l + fc+,3 + fc+) |fc +-^- fc °-° ' > 
with w = u>i + C1J2 — ^3 ; it is easy to check that this term is indeed non vanishing. Note also 
that A is the first non trivial contribution to the critical index rj of the Schwinger function 
of a d = 1 systems of interacting fermions. 

4.4 Schwinger functions. We will note repeat here the analysis of the Schwinger functions 
at the temperature scale, as one can proceed as in the d — 1 to obtain an expansion for the 
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Schwinger function once that the expansion for the effective potential is understood; see for 
instance [GM]. We only remark that the Aj and Ajj in (1.12) and (1.13) are indeed 0(X 2 ) 
as a consequence of 

J dk dk' + dk-g'}(ko,k' + + ujp Fl fc_) = J dk dk + dk'_g'}(ko, k+, k'_ + up F ) = (4.20) 



5. Conclusions 



5.1 Marginal Fermi liquids and Luttinger liquids. We can compare the behaviour of the 
half filled Hubbard model with cut-off with other models. We have found that the wave 
function renormalization has an anomalous flow up to exponentially small temperatures, 
Z g~ 1 = 1 + 0(A 2 ), sec (4.14); in the case of circular fcrmi surfaces one finds instead, see 
[DR], for |A| < e 

' = l + 0(e 2 7 *) (5.1) 



Z h 

which means that Zh = 1 + 0(A 2 ), up to exponentially small temperatures; the factor r ) h / 2 
in the r.h.s. of (5.1) is an improvement with respect to a power counting bound and is found 
by using a volume improvement based on the geometrical constraints to which the momenta 
close to the Fermi surface are subjected. An equation similar to (5.1) holds also for any 
symmetric smooth Fermi surfaces with non vanishing curvature; a proof can be obtained by 
combining the results of [BGM] with Appendix 2 of [DR]. 

The similarity of the equation for Zh with its analogous for one dimensional systems 
may suggest that the behaviour of the half-filled Hubbard model with cut-off up to zero 
temperature is similar to the one of a system of spinless interaction fermions in d = 1 (the 
so called Luttinger liquid behaviour). However this is false; in a Luttinger liquid in fact one 
has that 

Afc-i = X h + 0(e 2 ^) (5.2) 

a property known as vanishing of Beta function. One can easily check that this cancellation 
is not present in the half-filled Hubbard model with cut-off; in fact the dominant second 
order contribution to \h(ki,i, k-i,i, &!,/) containing only a = I internal lines is 



/ 



dk'f h (k , k+)f^ h {k Q , k + )-rJ-r S -H{t$ i Sin 2 fc_) 

' + 



[H(al sin 2 (fci _ + fc 2 ,_ - k-)X h X h - H{a 2 sin 2 (fc 3 _ - k 2 - + k-)X h X h ] (5.3) 

where the dependence from k of the A^ has not been explicitated. It is clear then that even 
at the second order the flow of Xh is quite complex, and we plan to analyze it in a future 
work, in order to understand the leading instabilities. Replacing H with 1 and having A^ 
not momentum dependent one recovers the d = 1 situation in which the beta function is 
vanishing. The theory resembles the theory of d = 1 Fermi systems in which each particle 
has an extra degree of freedom, the component of the momentum parallel to the flat Fcrmi 
surface, playing the role of a "continuous" spin index; and it is known that in d = 1 even a 
spin-1 index can destroy the Luttinger liquid behaviour. 

5.2 Marginal Fermi liquid behaviour close to half filling. A similar analysis can be per- 
formed in the case of the Hubbard model with cut-off close to half-filling {ji = — e with 
e small and positive); in such a case the Fermi surface is convex and with finite radius of 
curvature but still resembles a square with non flat sides and rounded corners. The propa- 
gator has the form .. ,„ x ^ z and it is easy to verify that, if 8 < C minf-1 where 

° — iko+2 cos cos /c_ — e J J ' ' L p J 

p is the radius of curvature of the Fermi surface, the bounds (3.3), (3.4) for the single scale 
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propagator still holds; the reason is that, up to temperatures greater than the inverse 
of the curvature radius, the bounds are insensitive to the fact that the sides of the Fermi 
surface are not perfectly fiat. One can repeat all the analysis of the preceding sections and 
it is found that the Schwinger functions behave like (1.12), (1.12) for A small enough and 
(3 < Cmin[min[l], e^ c l A ^ ]; in other words marginal Fermi liquid behaviour is still found 
close to half filling , up to such temperatures. 

On the other hand at lower temperatures, for min[l }</3< 1 (of course assuming 

nrin[i] < e( c l A l- 1 ) one can apply the results of [BGM] (valid for any convex symmetric and 
regular Fermi surface) so finding Z = 1 + C P [X 2 + 0(X 3 )] where C p is a constant which 
is very large for small e (and diverging at half filling e = 0). Hence, depending on the 
values of the parameters, one can have, in the low temperature region and before the critical 
temperature, two possibilities: the first is to have only marginal Fermi liquid behaviour 
Z = 1 + 0(A 2 log/3), and the second is to have marginal Fermi liquid behaviour up to 
temperatures 0(p~ 1 ) and then Fermi liquid behaviour up to the critical temperature. 



Acknowledgments This work was partly written during a visit at the Ecole Polytech- 
nique in Paris. I thank J. Magnen and V. Rivasseau for their warm ospitality and useful 
discussions. I thank also G. Benfatto and G. Gallavotti for important remarks. 



References 

[A] P.W.Anderson. The theory of superconductivity in high T c cuprates, Princeton 
University Press, Princeton (1997) 

[AGD] A.A.Abrikosov, L.P.Gorkov, I.Dzialoshinsky, Methods of quantum field theory in 
statistical physics, Prentice hall (1963) 

[BGM] G. Benfatto, A.Giuliani, V.Mastropietro. to appear Annales Poincare' 

[BM] G. Benfatto, V.Mastropietro. Renormalization group, hidden symmetries and ap- 
proximate Ward identities in the XYZ model. Rev. Math. Phys. 13 (2001), no. 
11, 1323-143 

[BoM] F.Bonetto, V.Mastropietro. Beta function and anomaly of the Fermi surface for a 
d = 1 system of interacting fermions in a periodic potential. Comm. Math. Phys. 
172 (1995), no. 1, 57-93. 

[DAD] S.Dusuel, F. De Abreu, B. Doucot Renormalization Group for 2D fermions with a 



flat Fermi surface cond-mat 0107548 



[DR] M.Disertori, V. Rivasseau. Interacting Fermi liquid in two dimensions. Comm. 

Math. Phys. 215, 251-290 and 291-341 (2000) 
[FMRT] J.Feldman, J. Magnen, V. Rivasseau, E. Trubowitz. An infinite volume expansion 

for Many Fermions Green functions Helvetica Physica Acta, 65, 679-721 (1992) 
[FSW] A.Ferraz, T.Saikawa, Z.Y. Weng Marginal fermi liquid with a two-dimensional 



Patched Fermi surface cond-mat9908111 
[FSL] J.O. Fjaerested, A. Sudbo, A. Luther. Correlation function for a two dimensional 

electron system with bosonic interactions and a square Fermi surface Phys. Rev. B. 

60, 19, 13361-13370 (1999) 
[GM] Gentile, .G, Mastropietro.V. Renormalization group for one- dimensional fermions. 

A review on mathematical results. Phys. Rep. 352 (2001), no. 4-6, 273-43 
[Le] A. Lesniewski: iJ/feciive action for the Yukawa 2 quantum field Theory. Comm. 

Math. Phys. 108, 437-467 (1987). 
[L] A.Luther, Interacting electrons on a square Fermi surface PRB, vol.30, n.16 (1994) 
[M] D.C.Mattis, Implications of infrared instability in a two dimensional electron gas 

PRB, Vol.36, n.l, 1987 
[Me] W.Metzner, Renormalization Groupanalysis of a two dimensional interacting elec- 
tron systemlnt. Jour. mod. Phys. 16,11, 1889-1898 (2001) 
[R] V. Rivasseau. The two dimensional Hubbard model at half-filling. I. Convergent 

contributions. J. Statist. Phys. 106 (2002), no. 3-4, 693-72 



1/ febbraio/ 2008; 23:12 



21 



[S] Z-X Shen et al, Science 267, 343 (1995) 

[Sh] R. Shankar, Renormalization Group approach to interacting fcrmions. Rev. Mod. 

Phys. 66 (1) 129-192 (1994) 
[VZS] CM. Varma, Z.Nussinov, W. van Saarloos Singular Fermi liquids cond-mat0103393] 
[VLSARp.M. Varma, P.B. Littlewood, S.Schmitt-Rink, E. Abrahams, A.E. Ruckestein Phe- 
nomenology of the normal state of the Cu-0 high T c superconductors Phys. Rev. 
Lett., 63, 1996, (1989) 

[VR] A.Viroszteck, J.Ruvalds. Nested fermi liquid theory. Phys Rev. B 42, 4064 (1990) 
[ZYD] A.T.Zheleznyak, V.M. Yakovenko, I.E. Dzyaloshinskii, Parquet solutions for a flat 
fermi surface Phys. Rev.B, 55, 3200, 1997. 



1/ febbraio/2008; 23:12 



22 



